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Abstract 

We study the cosmic expansion history by reconstructing the deceleration parameter q{z) 
from the SDSS-II type la supernova sample (SNIa) with two different light curve fits (MLCS2k2 
and SALT-II), the baryon acoustic oscillation (BAO) distance ratio, the cosmic microwave back- 
ground (CMB) shift parameter, and the lookback time-redshift (LT) from the age of old passive 
galaxies. Three parametrization forms for the equation of state of dark energy (CPL, JBP, and 
UIS) are considered. Our results show that, for the CPL and the UIS forms, MLCS2k2 SDSS-II 
SNIa+BAO+CMB and MLCS2k2 SDSS-II SNIa+BAO-FCMB-FLT favor a currently slowing- 
down cosmic acceleration, but this does not occur for all other cases, where an increasing cosmic 
acceleration is still favored. Thus, the reconstructed evolutionary behaviors of dark energy and 
the course of the cosmic acceleration are highly dependent both on the light curve fitting method 
for the SNIa and the parametrization form for the equation of state of dark energy. 
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I. INTRODUCTION 



One of the most important and mysterious issues of modern cosmology is the fact that 
our universe is undergoing an accelerating expansion . This observed phenomenon 
is usually attributed to the existence of an exotic energy component called dark energy 
(DE), which generates a repulsive force due to the negative pressure associated with it. 
The simplest candidate of DE is the cosmological constant with the equation of state 
(EOS) w = —1. If one generalizes the EOS w of DE from w = —1 to be an arbitrary 
constant w, the current astronomical observations js-S] show that w is confined to be 
|1 + wl < 0.06 at the la confidence level. However, the EOS might also be a func- 
tion of cosmic time. In order to unveil the evolutionary properties of dark energy using 
observations, one usually adopts a parametrization form with some free parameters for 
w{z), which may not be motivated by any particular foundamental theory. Examples of 
such kind are the Chevallier-Polarski-Linder (GPL) parametrization [8], the Jassal-Bagla- 
Padmanabhan (JBP) parametrization \^ as well as the Upadhye-Ishak-Steinhardt (UIS) 
parametrization 1^, and so on. Constraining the free parameters of a given parametriza- 
tion with observational data, one can obtain the evolutionary curve of w{z), which embod- 
ies the property of dark energy. For instance, some current observations give an indication 



a. 



that the EOS has crossed the phantom divider {w = —1) at least once 

Recently, Sahni et al. proposed a new diagnostic of DE, named Om{z) diagnostic. 
The advantage of this diagnostic, as opposed to the EOS, is that it depends on the 
first derivative of the luminosity distance diiz) fll|, and thus is less sensitive to the 
observational errors and the present matter energy density Qom- In addition, one can 
discriminate DE models with the EOS w < —l,w = —1 and w > —1 effectively by 
examining the slope of Om{z) even if the value of l^om is not exactly known, with positive, 
null, or negative slopes corresponding to w < —l,w = —1 or w > —1, respectively. 

Performing the Om(z) diagnostic {n] and analyzing the deceleration parameter q{z) 
with the Constitution type la supernova data (SNIa) [l2j and data from the baryon acous- 
tic oscillation (BAO) distance ratio of the distance measurements obtained at 2; = 0.2 and 
z = 0.35 in the galaxy power spectrum [3, Q by using the popular CPL parametriza- 
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tion, Shafieloo et al. found that Om{z) and q{z) increases markedly at the low redshifts 



z < 0.3 15j. This result suggests that the expansion acceleration of our universe is 
probably slowing down. However, the result obtained from a combination of the cosmic 
microwave background (CMB) shift parameter, SNIa and BAO is very well consistent 
with the ACDM model. So, there appears some tension between low redshift data (Con- 
stitution SNIa+BAO) and high redshift (CMB) one. Surprisingly, further analysis using 
a subsample (SNLS+ESSENCE+CfA) of the Constitution SNIa reveals that the outcome 
does not rely on whether the CMB data is added and the cosmic acceleration has been 
over the peak. It was therefore argued the tension could be either due to the systematics 
in some data or that the CPL parametrization is not versatile-enough to accommodate the 



evolution of DE implied by the data. This situation was also examined by Gong et al. [16 1 
recently through the reconstruction of Om{z). They found that both the systematics in 
data sets and the parametrization of DE influence the evolutional behavior of DE. It is 
worth noting that the results in Ref. [l5| are obtained only at the la confidence level and 
whether they are reliable at the 2a confidence level still needs to be checked. 

In this paper, we will re-examine this issue with three different parametrization forms 
for the EOS of DE (CPL, JBP and UIS js iO]), and 288 SNIa data points released by the 



Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey with two different light curve 
fits (MLCS2k2 fit and SALT-II fit). We hope this will help us further understand the 
influence of different light curve fitting methods ^. As in [l5|, the BAO distance ratio 
and the CMB shift parameter are also considered in our analysis. In addition, we use the 
lookback time-redshift (LT) from the age of old passive galaxies, since it has the advantage 
that the ages of distant objects are independent of each other, and thus may avoid biases 
that are present in techniques using distances of primary or secondary indicators in the 
cosmic distance ladder method. 



^ The influence of different light curve fits has been considered in jl6|. However, there, different data 
points are obtained with different light curve fits for Constitution SNIa. In the present paper, data 
points are the same for different light curve fits. 
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II. PARAMETRIZATION 



By assuming that the energy components of our universe are nonrelativistic matter 
and dark energy, the Friedmann equation can be expressed as 



DEi 



(1) 



where f^om is the current value of the dimensionless matter energy density, and VL^e is 
the dimensionless energy density parameter of dark energy, which can be expressed as 

l + wiz' 



^DE = (1 - ^Om) exp 



1 + z' 



-dz' 



(2) 



where w{z) is the equation of state of dark energy. 

Now, we consider three parametrization forms for the EOS. i.e., the GPL, JBP and 
UIS parametrization. The EOS w{z) for the GPL is js] 

WiZ 



w{z) = Wo + 



l + z 



Substituting it into Eqs. ([I1I2]), one has 
For the JBP parametrization, w{z) is js] 



exp 



31^12; 
l + z 



w{z) = Wq + 



WiZ 



[l + zf 



So we have 



E\z) = fiom(l + Zf + (1 - fio™)(l + Z^ 



|3(l+«;o) 



exp 



2(1 + z) 



(3) 



(4) 



(5) 



(6) 



And for the UIS parametrization 10|, the 'w{z) and E^{z) can be expressed respectively 
as, 

Wq + WiZ Z < 1, 
Wq + Wi Z >1, 



w{z) 



(7) 



and 



EHz) 



fiom(l + zf + (1 - l]o™)(l + 2)3(i+-o-»i) exp(3wi2) z < 1, 



fio™(l + zf + (1 - + ;2)3(l+-0+-i) 



(8) 



Z>1. 
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III. OBSERVATIONAL DATA 



The SNIa dataset used in our analysis is the 288 data points released by the Sloan 
Digital Sky Survey-II (SDSS-II) Supernova Survey l?!, which consist of 103 new SNIa 
from the first-year SDSS-II su per nova survey 17H19l|. 56_points from ESSENCE 201], 62 
from SNLS [sj, 34 from HST and 33 nearby SNIa [231. ^oi 288 SDSS-II SNIa data, 



two kinds of light curve fitting methods, i.e., the MLCS2k2 fit and the SALT-II fit, have 
been employed 17|. Here, for the sake of systematics influence check, we do our analysis 



of the SDSS-II data with both kinds of fits. 

The constraints from the SNIa data can be obtained by minimizing the following 
statistics: 



288 



xJ(/^o,p) = 



obs.i 



(9) 

i=l "A* 

where ^obs,i is the distance modulus estimated from the MLCS2K2 fit or SALT-II fit 
for the i'th supernova, Zi is its spectroscopically determined redshift, p stands for the 
complete set of model parameters, and fith{zi]p) is the theoretical distance modulus for 
a concrete cosmological model, which is calculated from 

Hthizi] /io, p) = 5 hg^Qidiizi] p)) + /io, (10) 



where represents the luminosity distance. For a flat universe, is given by 

1 



dL{z; p) 



1 + ^ 



dz'- , 
E{z\p) 



'111 



Let us note that here /io is a nuisance parameter. In order to marginalize over it, we 
expand (Eq. ©) with respect to /iq as 



xJ(/^o, P) = ^ - 2/ioS + /i^C 



where 



i?(p) 

C 



288 

E 

i=l 

288 

E 

i=l 

288 



err 



[/i 



obs', 



fithiziino = 0,p)] 



y-- 



(12) 
(13) 

(14) 
(15) 



Eq. ( !T2|) has a minimum at fiQ = B/C, and it is 

Xlip) = A{p) - ^ . (16) 



Thus, instead of minimizing xfiif^o,p), we can minimize xfi{p)> which is independent of 
/io, to obtain constraints on the model parameters. 

In the x^-expression, the distance-modulus uncertainty is given by 

= + «)' + «)^ (17) 

where a/** and a^"*(=0.16 for MLCS2k2 and 0.14 for SALT-II) are respectively the sta- 
tistical uncertainty and the additional (intrinsic) error derived from light-curve fitting 
method, o"™* is introduced to unitize P^r degree of freedom for the Hubble diagram 
constructed from the nearby SNIa sample and it is detailed in Ref. It has lit- 

tle effect on the results of model parameters, but improves the value of xt remarkably. 



For example, for SALT-II SNIa with and without cr^"*, the best fit values of the CPL 
parametrization are wq = —1.022 and wi = 0.090 with xfi = 268.66, and wq = —1.015 
and Wi = 0.200 with xf = 589.31, respectively, where fiom is given prior to be 0.262 



obtained from WMAP7 



271]. cr^ is the error which relates with the redshift uncertainty 



and can be computed from 

l + z 



InlOy ^(1 + ^/2)' ' ' 

where cx^ is the redshift uncertainty and is defined to be 

Here, az,spec and az^pec, as discussed in detail in Ref. jl?!], are the uncertainties from 
spectroscopic measurements and peculiar motion of the host galaxy, respectively. 

The BAO distance measurements used in our analysis are obtained at z = 0.2 and 



z = 0.35 from the joint analysis of the 2dF Galaxy Redshift Survey and SDSS data 14 1. 
The BAO distance ratio 

Dv{z = 0.35) 



Dviz = 0.20) 



1.736 ±0.065 (20) 
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(21) 



is a relatively model-independent quantity. Here Dv{z) is defined as 

For the BAO dataset, we can fit the model parameter p by performing statistics as 
follows 

2 . . _ \Dv{z = 0.35)/Dviz = 0.20) - 1.736]^ 
XbaoW - • ^^^) 

In addition to the low redshift SNIa and BAO, we add the high redshift CMB parameter 
which is the reduce distance at zis = 1090. The shift parameter 

Rip) = ^M^ - = 1.71 dz 0.019, (23) 



is used. We also apply 

0.0192 



- '^n^ (24) 



to the model parameter p for the CMB data. 

Besides the most common observational data sets above, we also perform an analysis 
combined with the lookback time-redshift data (LT), which is established by estimating 
the age of 32 old passive galaxies distributed over the redshift interval 0.11 < z < 1.84 and 
the total age of the universe f^^ [24]. The galaxy samples of passively evolving galaxies 
are selected with high-quality spectroscopy and the method used to determine ages of 
galaxy samples indicates that systematics are not a serious source of error for these high- 
redshift galaxies. In addition, this data set has the advantage that the ages of distant 
objects are independent of each other, and thus it may avoid biases that are present in 
techniques that use distances of primary or secondary indicators in the cosmic distance 
ladder method. As a result, these age data are different from the widely used distance 
one, such as SNIa, and it may help us gain more insight into the nature of dark energy. 
To estimate the best fit of model parameters, we minimize Xage 

T U jobs 



where, o"|i = af + a^^^, cTj is the uncertainty in the individual lookback time to the i*'* 
galaxy of the sample, CTiohs =0.7 Gyr stands for the uncertainty in the total expansion 



age of the universe (to^*), and r means the time from Big Bang to the formation of the 
object, which is the so-called delay factor or incubation time. Note that while the observed 
lookback time {t°^^{zi] r)) is directly dependent on r, its theoretical value {ti^Zi] p)) is not. 
Furthermore, in principle, it must be different for each object in the sample. Thus the 
delay factor becomes a "nuisance" parameter, we use the following method to marginalize 



over it 



25 



26| 



fOO / -j^ 

X^(P) = -21n / rfrexp ( - -x^ 



^age 



A- — + D-2hi 



erfcf , I 



2C \V2Cj 



(26) 



where 



A 



A 2 A 1 

D is the second term of the rhs of Eq. fl25l) , 

^ = tL{z,-^)-[tt -t{zdl (28) 
and erfc(x) is the complementary error function of the variable x. 

IV. RESULTS 

We first explore three popular parametrization forms by using SDSS-II SNIa, BAO 
and CMB. A comparison of two light curve fits of SDSS-II SNIa is given. The results are 
shown in Figs. ([11 [2]). 

Fig. ([T]) gives the contour diagrams on wq—Wi panel at the 68.3% and 95.4% confidence 
levels for three parametrization forms. Since different data sets give different best fit 
values of Qom; set prior i^om = 0.262 in this paper, which is the best fit value from the 



WMAP7 271]. The top, middle, and bottom panels are the results of the CPL, JBP and 
UIS, respectively. The big red points denote the spatially flat ACDM model {wq = — 1 
and Wi = 0). In the left panels, the SNIa is obtained with the SALT-II light curve 
fit, whereas, in the right panels, it is given with the MLCS2k2 one. The blue dotted, 
yellow dot-dashed, green dashed and pink solid lines show the results from SDSS-II SNIa, 
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BAO+CMB, SDSS-II SNIa+BAO, and SDSS-II SNIa+BAO+CMB, respectively. From 
this figure, it is interesting to see that, for all three different parametrization forms, the 
ACDM model is consistent with the SNIa, SNIa+BAO and SNIa+BAO+CMB at the la 
confidence level when the SALT-II fit is used. This consistency is, however, broken at 
the 2cr confidence level when the fit is changed to MLCS2k2. Thus, the consistency of 
the ACDM model with data depends crucially on the type of fit used. This agrees with 



what was obtained in Ref. 



from the SDSS-II SNIa with the CPL parametrization. 



Let us note that this discrepancy between two analysis methods has been pointed out in 



Refs. 



paper \v\. Furthermore, we find that the SDSS-II SNIa with SALT-II fit is consistent 
with BAO+CMB very well. However, when the MLCS2k2 fit is used, there exists a 
tension between SNIa and CMB+BAO at the la confidence level. Actually, this tension 
also exists between other SNIa sets. Gold, for an example, and CMB+BAO 31|- 

Fig. ([2]) shows the evolutionary curves of the deceleration parameter q{z) recon- 
structed from the SDSS-II SNIa (SALT-II and MLCS2k2)+BAO and the SDSS-II 
SNIa+BAO+CMB for three parametrization forms. The gray regions and the regions 
between the two dashed lines represent the Icr confidence level for q{z) obtained from 
SNIa+BAO+CMB and SNIa+BAO, respectively. The dot-dashed lines show the best fit 
curves of the spatially fiat ACDM model. The left panels in Fig. (j2]) reveal that, when the 
SALT-II fit is considered, independent of the parametrization forms, both SNIa+BAO and 
SNIa+BAO+CMB support an accelerating cosmic expansion and the acceleration seems 



and it was also found with a simple wCDM model in the original SDSS-II 
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where it was 



to be speeding up. This is different from the results obtained in Refs. 
found that the Constitution SNIa+BAO favor a slowing-down of the cosmic acceleration 
at the low redshifts. However, the right panels show that, for the MLCS2k2 fit, the best- 
fit results obtained from SNIa+BAO+CMB for the CPL and UIS parametrization forms 
favor a slowing-down of the cosmic acceleration. While, when the JBP parametrization 
form is used, the results obtained from SNIa+BAO+CMB change markedly and an in- 
creasing cosmic acceleration is favored. If the CMB data is not included, the results are 
consistent for three parametrization forms and observations favor an increase of the cos- 
mic acceleration. In addition, we also find that the SALT-II SDSS-II SNIa gives a best-fit 
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current acceleration which is larger than that from the ACDM, whereas the MLCS2k2 
SDSS-II SNIa gives a one which is less. Thus, different light curve fits of SNIa may yield 
completely different behavior of dark energy and the parametrization forms also matter. 
For the purpose of unfolding the uncertainty in supernova light curves, we plot the differ- 
ence of the distance modulus obtained with two fits for each supernova in Fig ([3]). From 
this figure, it is easy to see that for most low redshift supernova data points {z is about 
less than 0.5) their distance moduli from the SALT-11 fit are smaller than those from 
the MLCS2k2 fit, which apparently lead to different results for the cosmic evolution from 
SNIa. 

Since the LT has the advantage of avoiding the biases existing in data obtained from the 
cosmic distance ladder method, we add it in our discussion to obtain more information on 
the evolutional behavior of dark energy. The results are shown in Figs. (jH El [6]). Fig (jl]) 
gives the allowed region of the model parameters Wq and Wi at the 68.3% and 95.4% 
confidence levels. The dotted, dashed and solid lines represent the results from SNIa, 
BAO+CMB+LT and SNIa+BAO+CMB+LT, respectively. Comparing the right panels 
of Fig. ([T]) and Fig. (jlj), one can see that when LT data is added, the tension between 
the MLCS2k2 SDSS-II SNIa and other data sets becomes more severe and it exists at the 
2(7 confidence level. However, LT data render the results to be more consistent with the 
ACDM, since MLCS2k2 SDSS-II SNIa+BAO+CMB+LT allow the ACDM, but MLCS2k2 
SDSS-II SNIa+BAO+CMB rule out it at the 2a confidence level. 

Fig. ([5]) gives the evolutionary curve of q{z) reconstructed from 
SNIa+BAO+CMB+LT. The gray regions show 68.3% and 95.4% confidence levels, 
the thick solid lines are the best fit results, and the dot-dashed lines indicate the 
spatially fiat ACDM. We find that the results are similar to that obtained from SDSS-II 
SNIa+BAO+CMB. Only in the cases of the CPL and UIS parametrizations and the 
SDSS-II SNIa with MLCS2k2 fit, do SNIa+BAO+CMB+LT favor that the cosmic 
acceleration is slowing down. In all other cases, an increasing of the cosmic acceleration 
is favored. That is to say, the inclusion of the LT data does not change the evolutionary 
behavior of q{z) markedly. 

To further confirm quantitatively the behavior of q{z), we plot, in Fig. (Q, the 
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evolutionary curves of the jerk parameter j{z), which is proportional to the deriva- 
tive of q{z). From the left panel of this figure, one can see that SALT-11 SDSS-II 
SNIa+BAO+CMB+LT favor an increasing cosmic acceleration at the present since the 
value of j(0) is largely in the region of > 0, but they cannot rule out the possibility 
of a decrease of the cosmic acceleration since j(0) < is still allowed at the 2a confi- 
dence level. The right panel of Fig. (E]) shows that although, for the CPL and the UIS 
parametrizations, observations favor a slowing down of the cosmic acceleration, they still 
allow a currently increasing cosmic acceleration. For the JBP parametrization, similar re- 
sults from SDSS-II SNIa with SALT-II and MLCS2k2 fits are obtained, and observations 
always favor an increasing cosmic acceleration. 



V. CONCLUSION 



We focus, in this paper, on analyzing the property of dark energy by reconstructing 
the deceleration parameter q{z) from the SDSS-II SNIa with two different light curve fits 
(SALT-II and MLCS2k2) along with BAO, CMB and LT. Three different parametrization 
forms of the EOS of dark energy (CPL, JBP and UIS) are investigated. We find that, 
when the SALT-II light curve fit is considered, independent of the parametrization forms 
SNIa+BAO, SNIa+BAO+CMB as well as SNIa+BAO+CMB+LT favor an accelerating 
cosmic expansion which is currently speeding up. We also find that there is an tension 



T). T 
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lese differ from 
where it was 



between SNIa and other data sets (CMB+BAO or CMB+BAO+ 
what was obtained from the Constitution SNIa+BAO in Refs. 
found that the cosmic acceleration is probably slowing down now. However, when the 
MLCS2k2 light curve fit is considered, the results are dependent on the parametrization 
forms. For the CPL and UIS parametrization forms, the best fit results obtained from 
SNIa+BAO+CMB and SNIa+BAO+CMB+LT show that a slowing-down acceleration of 
the cosmic expansion is favored, and this, however, does not occur for SNIa+BAO. For 
the JBP parametrization, an increasing cosmic acceleration is always favored. Finally, 
from the evolutionary curves of j{z), we find that both an increasing and a decreasing 
of the current cosmic acceleration are allowed by the SNIa+BAO+CMB+LT at the 2a 
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confidence level. 

Thus, the evolutional behavior of dark energy reconstructed and the issue of whether 
the cosmic acceleration is slowing down or even speeding up is highly dependent upon the 
SNIa data sets, the light curve fitting method of the SNIa, and the parametrization forms 

of the equation of state. In order to have a definite answer, we must wait for data with 
more precision and search for the more reliable and efficient methods to analyze these 
data. 
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FIG. 1: The 68.3% and 95% confidence level regions for wq versus wi. The blue dotted, yel- 
low dot-dashed, green dashed and pink solid lines represent the results obtained from SNIa, 
BAO+CMB, SNIa+BAO and SNIa+BAO+CMB, respectively. The left panels show the re- 
sults obtained from SNIa with SALT fit, while the right panels are the results from SNIa with 
MLCS2k2 fit. The red point at z(;o = — 1, 't^i = represents the spatially flat ACDM model. 
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FIG. 2: The gray regions and the regions between two dashed hnes show the evolutionary be- 
haviors of q{z) at the 68.3% confidence level obtained from SNIa+BAO+CMB and SNIa+BAO, 
respectively. The left panels show the results obtained from SNIa with SALT fit, while the right 
panels are the results from SNIa with MLCS2k2 fit. The dot-dashed lines represent the best-fit 
spatially flat ACDM model. 
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FIG. 3: The difference in the distance modulus obtained with two fits for each supernova. The 
subscripts s and m represent the SALT-II fit and the MLCS2k2 fit respectively. 
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FIG. 4: The 68.3% and 95% confidence level regions for wq versus wi. The blue dotted, green 
dashed and pink solid lines represent the results obtained from SNIa, BAO+CMB+LT and 
SNIa+BAO+CMB+LT, respectively. The left panels show the results obtained from SNIa with 
SALT fit, while the right panels are the results from SNIa with MLCS2k2 fit. The red point at 
wq = —1, wi = represents the spatially fiat ACDM model. 
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FIG. 5: The gray regions show the evolutionary behaviors of at the 68.3% and 95% confi- 
dence levels obtained from SNIa+BAO+CMB+LT. The dot-dashed lines represent the best-fit 
spatially flat ACDM model. The left panels show the results obtained from SNIa with SALT 
fit, while the right panels are the results from SNIa with MLCS2k2 flt. 
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FIG. 6: The gray regions show the evolutionary behaviors of j(z) at the 68.3% and 95% con- 
fidence levels obtained from SNIa+BAO+CMB+LT. The left panels show the results obtained 
from SNIa with SALT fit, while the right panels are the results from SNIa with MLCS2k2 fit. 
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